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ABSTRACT 

We present results of a survey of Lyman a emitters (LAEs) at z = 6.5 which is 
thought to be the final epoch of the cosmic reionization. In a k 530 arcmin 2 deep 
image of the SSA22 field taken through a narrowband filter NB912 installed in the 
Subaru/Suprime-Cam, we have found only 14 LAE candidates with LLya > 3 x 10 42 
erg s . Even applying the same colour selection criteria, the number density of the 
LAE candidates is a factor of 3 smaller than that found at the same redshift in 
the Subaru Deep field (SDF). Assuming the number density in the SDF is a cosmic 
average, the probability to have a number density equal to or smaller than that found 
in the SSA22 field is only 7% if we consider fluctuation by the large-scale structure 
(i.e. cosmic variance) and Poisson error. Therefore, the SSA22 field may be a rare 
void at z — 6.5. On the other hand, we have found that the number density of i'-drop 
galaxies with 25.5 < z' < 26.0 in the SSA22 field agrees well with that in the SDF. 
If we consider a scenario that a larger neutral fraction of intergalactic hydrogen, cchi, 
in the SSA22 field obscures a part of Lya emission, xm in the SSA22 field should be 
about 2 times larger than that in the SDF. This can be translated into xm < 0.9 at 
z = 6.5 in the SSA22 field. A much larger survey area than previous ones is required 
to overcome a large fluctuation reported here and to obtain a robust constraint on xm 
at the end of the reionization from LAEs. 
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1 INTRODUCTION 

Survey of Lya emitters (LAEs) at z > 6 has been 
significantly developing i n the last decade since the 
initial discovery papers iHu. McMahon fc Cowiel Il999l; 



initial discovery papers (|llu, iviclvianon, & Uowic lyyy: 
Rhoads fc Malhotral l200ll : iKodaira et al.l 120031 ; I Ive et al.l 



2006) 



At 



5.7, many LAE surveys wit h the 
Subaru/Suprimc-Cam (S-Cam; Mivazaki ct al. 2002) have 



been made so far JHu et alj |2004 lAiiki et ajTl200l l200d: 



Shimasaku et alj [20061: iMuravama et alj 120071 ; lOuchi et al.l 
20081 ). In particular. Shimasaku et al. ( 20061 ) have performed 



a deep LAE survey at z = 5.7 in the Subaru Deep Field 
(SDF) and present a semina l luminosity function (L F) of 
Lya emission at the redshift. IMuravama et all (J2007) have 
performed a very wide (~ 2 deg 2 ) survey in the COS- 
MOS field and found no evidence of a large-scale cluster- 



* This work is based on data collected at Subaru Telescope and 
obtained from the SMOKA, which is operated by the Astronomy 
Data Center, National Astronomical Observatory of Japan, 
f E-mail: akinoue@las.osaka-sandai.ac.jp 



ing of brigh t (L Lva > l x 10 43 erg s" 1 ) LAEs at z = 
5.7. Finally, lOuchi et alj (|2Q08T > have reported results of a 
wide (« 1 deg 2 ) and deep survey in the Subaru/XMM- 
Ne wton Deep Survey (SX DS) field an d confirmed the LF 
bv lShimasaku et all (|2006h . In addition. lOuchi et al.1 (|2008l ) 
report a factor of 5 variation of the number density of LAEs 
among their 5 fields of view of S-Cam. 

On the other h and, few surveys had been reported 
to date at z ~ 6 5 jTaniguchi et al.l 120051 : iHu et al.ll2005l : 
Hu & Cowic 2006) an d a significant improvement is aris- 
ing (|Ouchi et al.ll2010l : IHu et ; alj|2010h . A largest sample of 
LAEs at z = 6.5 reported so far is the one in the SDF 
presented by [Taniguchi et al.l |2005l ). Based on this sample, 
Kashikaw a et al.l (|2006h have constructed the Lya LF at the 
redshift and found that a significa nt reduction at the bright - 
end relative to the z = 5.7 LF bv lShimasaku et al.1 (2006). 
They interpret this reduction as a signature of the end epoch 
of the cosmic reionization. 



The cosmic reionization is one of the most drastic event 
in the history of the universe and it happened at z ~ 10 
jPunklev et al.1 [2009 ) . The reionization was probably trig- 
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gered by the first-generation of stars and gaiaxies. Thus, 
the event contains the information of the structure forma- 
tion in the early universe. The subsequent galaxy formation 
proceeded in the reionized universe and was significantly af- 
fected by the ionizing background radiation. This is a com- 
plex astrophysical process (or feedback) which is one of the 
most import ant issues to be resolve d in observational cos- 
mology (e.g. jLoeb fc Barkanall200lh . 

The discovery o f the Gunn-Peterson trough 
l)Gunn fc Peterson! 1 19651) in spectra of QSOs a nd a GRB 
at z ~ 6 | Becker et al.ll200ll ; pTotani et alJbOOfl) marks the 
end o f the reionizatio n epoch (e.g., iFan. Cari lli. fc Keating! 
2006i). Furthermore, IKashikawa et alJ \200& ) argue that 
the evolution of Lya LF between z — 5.7 and z — 6.5 
which they found is another signature of the end epoch; 
since the measurements of ultraviolet (UV) LF of the LAEs 
agree very well between the two redshifts, the apparent 
faintness of Lya emission of z = 6.5 LAEs is caused by an 
increase of attenuation by neutral hydrogen remained in 
the intergalactic medium (IGM) at z > 6. 

A significant variation in the IGM transmission found 
in the sp ectra of z ~ 6 QSOs suggests inhomogeneous reion- 
ization {Diorgovski. Bogosavlievic. fc Mahaball [2006) . Nu- 
merical simulations a lso show inhomogeneity of the reion- 
ization process (e.g.. IcTardi, Stoehr. fc White! [2003) . This 
means that there may be a large variation in the end epoch 
of the reionization among lines of sight. If it is true, the 
number density of z = 6.5 LAEs may variate signif i cantly 
among different survey fields. Indeed, iHu fc Cowiei (|2006T ) 
report such a large variation between two fields of view of 
S- Cam. Therefore, anoth er survey of z = 6.5 LAEs similar 
to lTaniguchi et all (|2005l ) is indispensable to discuss such a 
field-to-field variation. 

This paper presents a photometric sample of z = 6.5 
LAEs found in the archival imaging data taken with S-Cam 
towards the SSA22 field. Some of th e LAE s have been con- 
firmed by spectroscopy in IHu et all (|2010l ). Then, we com- 
pare the number density of the LAEs in the SSA22 field 
with that in the SDF in order to discuss the field-to-field 
variation quantitatively. The structure of the rest of this pa- 
per is as follows; first, we describe the imaging data used in 
this paper in §2. Then, we present the photometric sample 
of z = 6.5 LAEs in the SSA22 in §3. In §4, we compare the 
Lya LF derived in the SSA22 with that in the SDF, and 
also compare the number density of i'-drop galaxies in the 
two fields. In §5, we discuss the cause of the small num- 
ber of LAEs in the SSA22. The final section provides the 
conclusions of this paper. 

The cosmology adopted in this paper is a flat ACDM 
cosmology with h — 0.7, Qm = 0.3, and Qa = 0.7. All mag- 
nitudes are described by the AB system: AB = — 2.5 log — 
48.60, where f v in unit of erg s _1 cm - 2 Hz" 1 (|OkeHl974T ). 



2 IMAGING DATA 

We have collected imaging data of the SSA22 field 
(a = 22 h 17 m 34 s , 5 = +00° 1 7' 00"; J2000) taken 
with Subaru/Suprime-Cam (S-cam; iMivazaki etafl l2002h 
from the archive data system, SMOKA (Subar u-Mitaka- 
Okayama-Kiso Archive System; iBaba et al.ll2002T ). The col- 
lected broadband data are B, V, R, i' , and z'. In order to 
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Figure 1. Efficiency curves of the narrowband filters to select 
z ~ 6.5 LAEs. NB912 was used for the SSA22 field (this work) 
but NB921 was used for the SDF (Taniguchi et al. 2005). The 
efficiencies include the CCD quantum efficiency as well as the 
filter transparencies. 



select z ~ 6.5 LAEs, we use the imaging data through a 
narrowband NB912 (central wavelength of 9,139 A and full 
width at half maximum [FWHM] of 134 A) whose efficiency 
curve is shown in Figure 1. Table 1 is a summary of the 
imaging data. 

All the data were reduced using sdfred (Yag i et"aD 
l2002l ; lQuchi 7 et al . 2004) and iraf with the standard manner. 
In this paper, we selected only frames taken under a good 
seeing condition. Then, FWHMs of stellar objects in the fi- 
nal images are relatively good as summarised in Table 1, 
except for R band. The relative accuracy in astrometry in 
the final images is less than 0.5 pix (pixel scale is 0."202). 
The flux density calibrations of z and NB912 images were 
done by a photometric standard star GD 50 and a spectro- 
scopic standard star G 93-48, respectively. The calibrations 
of other broadband images were done so as that the colour 
distributions of z around 25.5 AB objec ts in the SSA22 
agree with those in the SDF reported by IKashikawa et all 
(2004), af ter the correction of the Galactic d ust extinction 
based on lSchleeel. Finkbeiner. fc Davisl (|l998l ). 

The data of z' and NB912 were taken with two posi- 
tion angles with a difference of 90°. While the S-cam has 
a field of view of 34' x 27', therefore, the area with a good 
homogeneous quality in the final co-added image is reduced 
to 22. '7 x 23. '6 for the two filters. We use only this area to 
search for LAEs at z ~ 6.5. The limiting magnitudes of z 
and NB912 in Table 1 are also measured in this area. The 
co-moving depth of our survey is 41.6 hg \ Mpc since NB912 
effectively captures Lya emission line between z — 6.46 and 
z = 6.57. Finally, our survey volume is 1.3 x 10 5 Mpc 3 
in co-moving unit. 

As a comparison field, we adopt the SDF (a = 
13 fe 24 m 39 s , 8 = +27°29'25"; J2000). IKashikawa et all 
(2004) published the reduced and calibrated images of B, 
V, R, i', z', and two narrowband NB816 and NB921. The 
NB921 has the central wavelength of 9,196 A and FWHM of 
the transmission of 132 A, whi ch are quite similar to those 
of NB912 as shown in Figure l. lTaniguchi et al.l {2005) and 
IKashikawa et all (|2006l ) reported their results of LAE survey 
using NB921 in the SDF. We use these imag e s and photo- 
metric catalogue reported in IKashikawa et"aL (|2004l ) and in 
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Table 1. Imaging data of the SSA22 field used in this paper. 



Filter 


Observed year (P.I.) 


Exposure time (h) 


PSF FWHM (") 


1-cr limit (AB) a 


B 


2008 (Y. Nakamura) 


2.6 


0.78 


29.08 


V 


2002 (T. Hayashino), 2003 (P. Capak), 2008 (Y. Nakamura) 


2.1 


0.82 


28.85 


R 


2000, 2001 (E. Hu) 


2.9 


1.06 


28.22 


i' 


2001 (P. Capak) 


1.4 


0.78 


28.28 


z' 


2001 (Y. Komiyama), 2002, 2003 (P. Capak) 


2.9 


0.77 


27.77 


NB912 


2002, 2003 (P. Capak) 


12.7 


0.77 


27.79 



a Aperture diameter is twice of the PSF FWHM. The Galactic dust extinction is corrected based on Schlegel et al. (1998). 
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Figure 2. z' number counts in the SSA22 field (plus) and SDF 
(cross). 

iTaniguchi et al. (2005) for a comparison with our LAE can- 
didates in the SSA22 field in §4. 

Before selecting LAEs, we should confirm the accuracy 
of the flux calibration of our data in the SSA22 field. For 
this aim, we compare the z' number count in the SSA22 
field with that in the SDF. The result is shown in Figure 2. 
We have found an excellent agreement in 22 < z' < 26. The 
differences in brighter and fainter magnitudes are due to 
saturation of counts and detection limits. We have also con- 
firmed an excellent agreement in NB number counts. This 
ensures the accuracy of our flux calibration in the SSA22 
field. 



3 SELECTION OF Z — 6.5 LAE CANDIDATES 

First, we mad e object detection in th e final NB912 image 
by SExtractor (|Bertin fc Arno uts 1996) with the criterion of 
"5 connected pixels above 2-a" . We detected 64,787 objects. 
Then, we selected objects satisfying all the following criteria: 

(a) 22.00 < NB912 < 26.04 (5cr), 

(b) z' - NB912 > 1.0, 

(c) z> - NB912 > -2.51og 10 (l - 3a[F z , /F NBgi2 ]), 

(d) »'-*'> 1.3 or %' > 27.53 (2a), 

(e) B > 27.89, V > 27.66, and R > 27.03 (3er), 

where a indicates 1-a uncertainty of photometry in each 
band and cr[i ? 2 //J 7 NB9i2] means 1-cr uncertainty of the flux 
density ratio of z' and NB912. The criterion (c) ensures 
that the colour excess deviates from a flat F v spectrum by 
more than 3-cr of the flux density ratio. These criteria are 
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Figure 3. Colour-magnitude diagram between z' — NB912 and 
NB912 in the SSA22 field. All objects detected in NB912 arc 
shown by small dots. The dot-dashed and dashed lines show the 
applied colour and magnitude cuts and 3-cr uncertainty in the flux 
density ratio of z' and NB912, respectively. The selected LAE 
candidates are shown by filled squares. The open squares are one 
additional LAE candidate and one spectroscopic LAE whose pho- 
tometry is affected by a neighbouring object. If objects are not 
detected in z' , the flux densities are replaced by 1-cr uncertainty 
in the band. 

equiva lent with those in the SDF applied bv lTaniguchi et all 
(2005), but our B, V, R, and i depths are 0.2-0.8 AB shal- 
lower than those in the SDF. The resulting number of ob- 
jects is 12, after visual check to remove objects evidently 
affected by spiders of bright foreground stars. We show 
the colour-magnitude diagram, the spatial distribution, and 
thumbnail images of the LAE candidates in Figures 3, 4, 
and 5, respectively. The coordinate and photometric mea- 
surements of these objects are summarised in Table 2. 

In addition, we found one object which may be a LAE 
at z ~ 6.5 in a careful visual checking of hundreds objects 
satisfying the magnitude and colour cuts (criteria [a]-[c]) 
but not the drop-out criteria ([d] and [e]). The photometry 
in shorter wavelengths of this object is probably affected by 
a bright neighbour as found in Figure 5 second line from 
the bottom. If the effect by the neighbour was removed, this 
object would satisfy all the selection criteria. As shown in 
Figure 3 by an open square (brighter one), this object has 
an excess large enough in z' — NB912. We add this object as 
No. 13 in the final list of our LAE candidates. 

Furthermore, we add another o bject as No. 14 i n the 
LAE list. This object is reported in iHu et all |2010) as a 
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Table 2. Observed properties of z ~ 6.5 LAE candidates in the SSA22 field. 



Object 




RA 


Dec 


i' a 


z' a 


NB912 a 


Remarks 








(J2000) 


(J2000) 


(AB) 


(AB) 


(AB) 






SSA22- 


z6p5-l 


22:17:42.84 


+00:18:08.5 


> 28.28 


26.57 


25.16 


HC221742-001808 (z ap 


= 6.4692) 


SSA22- 


z6p5-2 


22:16:44.62 


+00:05:54.8 


27.77 


> 27.77 


25.39 






SSA22- 


z6p5-3 


22:18:10.79 


+00:04:02.0 


> 28.28 


> 27.77 


25.51 






SSA22- 


z6p5-4 


22:17:40.13 


+00:08:14.1 


27.61 


27.22 


25.74 






SSA22- 


z6p5-5 


22:16:45.15 


+00:12:38.3 


> 28.28 


27.72 


25.76 






SSA22- 


z6p5-6 


22:16:51.72 


+00:04:53.5 


27.60 


27.32 


25.77 






SSA22- 


z6p5-7 


22:18:07.05 


+00:07:26.9 


> 28.28 


> 27.77 


25.81 






SSA22- 


z6p5-8 


22:17:57.12 


+00:18:19.4 


> 28.28 


27.21 


25.81 






SSA22- 


z6p5-9 


22:17:22.73 


+00:21:44.2 


28.23 


> 27.77 


25.84 






SSA22- 


z6p5-10 


22:17:23.69 


+00:10:23.0 


> 28.28 


27.55 


25.84 






SSA22- 


z6p5-ll 


22:17:38.28 


+00:09:09.1 


28.23 


> 27.77 


25.85 


HC221738-000909 (z sp 


= 6.4811) 


SSA22- 


z6p5-12 


22:18:01.65 


+00:22:21.1 


> 28.28 


> 27.77 


25.91 


HC221801-002220 (z sp 


= 6.5360) 


SSA22- 


z6p5-13 b 


22:18:01.41 


+00:14:12.2 


27.15 


25.99 


24.71 






SSA22- 


z6p5-14 b 


22:17:25.65 


+00:11:18.4 


27.16 


26.19 


24.74 


HC221725-001119 (z sp 


= 6.5142) 



a Magnitudes measured within a circular aperture diameter of twice of PSF FWHM of each band. Magnitudes below the 1-<t limit are 
denoted by the lower limit. 

b Photometry of this object is affected by a neighbouring object. 
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Figure 4. Spatial distribution of candidates of galaxies at z > 6 
in the SSA22 field. The filled squares are z ~ 6.5 LAE candidates 
and the open squares are one additional candidate and one spec- 
troscopic LAE whose photometry is affected by a neighbouring 
object. The circles are i'-drop galaxies with 25.5 < z' < 26.0. The 
dashed line show the survey area. 



spectroscopic LAE. As shown in Figure 5 bottom line, this 
object has a close neighbour which is probably low-z ob- 
ject. The separation is very close (~ 1"), so that we missed 
finding it. The colour excess z'— NB912 of this object satis- 
fies with the criterion (b) shown in Figure 3 (fainter open 
square) . 

In our photometric LAE sampl e, there are 4 o bjects 
confirmed to be LAEs at z = 6.5 bv lHu et al.l (|2010l ) spec- 
troscopically. We note their redshift in Table 2. Other ob- 
jects are possibly foreground objec t s like [O n], [O ill], or Ha 
emitters. Indeed, frvashikawa et al.l (|2006l ) have found one [O 
ill] emitter and five unidentified single emission line objects 
among 22 LAE candidates which they took the spectra. The 
five single emission line objects are either LAE or [O il] 
emitter. The contamination fraction in our LAE candidates 
would be similar to that in the SDF sample because the se- 



lection criteria are essentially the same. However, this issue 
should be confirmed by a follow-up spectroscopy in future. 



4 COMPARISON WITH SUBARU DEEP 
FIELD 

We have found 14 candidates of LAE at z ~ 6.5 in the 
SSA22 field. In this section, we compare statistics of our 
LAE c andidates with those in an other large survey in th e 
SDF bv lTaniguchi et all i|2005l ) and lKashikawa et all (2006). 



4.1 Lya luminosity function 

First, we compare the cumulative LFs of the LAE candi- 
dates in the SSA22 field and in th e SDF. While we found 
only 14 candidates in the SSA22, [Taniguchi et all (|2005f ) 
reported 57 candidates in the SDpJ Among their 57 ob- 
jects, one object was found to be a foreground [O ill] emit- 
ter by follow-up spectroscopy and other five objects may 
be also foreground co ntaminations l|Taniguchi et all 120051 ; 
iKashikawa et al.ll2006l) . However, we compare our 14 can- 
didates with their 57 candidates because we may also have 
foreground contamination in our 14 photometric samples. 

We estimat e Lyq flux by the same method adopted in 
iTaniguchi et~ail l|2005l ): 



rcont 



AAiv 



2/2, 



(1) 



/ LyQ = /r s (NB912)AA NB 9i2 

where /£ bs (NB912) is the observed NB912 flux density (per 
A), /^ ont is the continuum flux density, and AAnb912 = 134 
A is FWHM of the NB912 efficiency curve. The division 
by 2 accounts for the continuum break at the shortward of 
Lya emission line because of intergalactic attenuation. The 
continuum flux density is estimated from 



1 They published 58 objects but one object, No. 22, has z' — 
NB921 = 0.87 which does not satisfy their criterion z 1 - NB921 > 
1.0. Thus, we remove this object from discussions in this paper. 
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Table 3. Estimated properties of z : 
SSA22 field. 
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Figure 5. Thumbnail images of the z ~ 6.5 LAE candidates. 
North is up and east is left. The field of view is 9" X 9". 



f cont = {f o hs{z , )AKi _ / obB (NB912)AANB912}/AA cff ( (2) 

where f\ bs (z') is the observed z' flux density, AA Z / = 1180 
A is the band-width of the %' filter, and AA°? = 450 A 
is an effective band-width for the continuum, which is the 
width between the longward edges of z and NB912 band- 
widths. Note that we o mit a numerical factor m ultiplied to 
/£ bs (NB912)AA NB 9i2 in iTaniguchi et all |2005) because it 
is actually the efficiency of the z' filter at Lya wavelength 



Object 






EW obs c 




HO 42 pro- Q~ ^\ 
^1U erg b ) 


HO 28 pro- s — 1 H7~ x 1 

l±W Cl g O .1.1 A. J 


cm 2 


SSA22-z6p5-l 


5.74 


8.5 


2.4 


SSA22-z6p5-2 


5.91 


< 4.8 


> 4.5 


SSA22-z6p5-3 


5.26 


< 4.8 


> 4.0 


SSA22-z6p5-4 


4.25 


< 4.8 


> 3.2 


SSA22-z6p5-5 


4.20 


< 4.8 


> 3.2 


SSA22-z6p5-6 


4.15 


< 4.8 


> 3.1 


SSA22-z6p5-7 


3.99 


< 4.8 


> 3.0 


SSA22-z6p5-8 


3.09 


4.8 


2.3 


SSA22-z6p5-9 


3.89 


< 4.8 


> 2.9 


SSA22-z6p5-10 


3.89 


< 4.8 


> 2.9 


SSA22-z6p5-ll 


3.85 


< 4.8 


> 2.9 


SSA22-z6p5-12 


3.66 


< 4.8 


> 2.8 


SSA22-z6p5-13 d 


8.16 


< 16 


> 1.9 


SSA22-z6p5-14 d 


8.49 


< 12 


> 2.6 



a Lyct luminosity estimated by equation (1). 

b UV continuum luminosity estimated by equation (2). 

c Observed equivalent width of Lya. 

d The derived UV luminosities are an upper limit because z' flux 
densities are contaminated by a neighbouring object. 



Table 4. The observed number density in unit co-moving volume 
of z ~ 6.5 LAE candidates in SSA22 and SDF. 



L hya Af(SSA22) a N(SDF) a 

(10 42 erg s" 1 ) (10 -5 Mpc" 3 ) (10~ 5 Mpc~ 3 ) 



> 10 


Q+1.4 


o 0+1.6 

z -°-i.o 


> 7 


i -°-1.0 


s.i±?:g 


> 5 


n 7+2.5 
d - '-1.6 


12.9l 2 j 


> 3 


10.4+»;» 


25 - 8 ±i:S 



a 1-cr uncertainties are based on Poisson statistics (Gehrels 1986). 

relative to the average efficiency of the filter and about unity. 
Although this modification is not very important, we re- 
calculate Lya fluxes of the LAE candidates in the SDF. 
When the estimated /^ ont is less than its 1-cr uncertainty 
estimated from an error propagation in equation (2), we set 
jcont _ o m e q U ation (1). Finally, we converted /Lya into 
Lya luminosity which are summarised in Table 3. In the 
table, we also list the continuum luminosity density and the 
observed equivalent width of Lya for which we adopted the 
estimated 1-cr uncertainty as an upper or lower limit. 

Figure 6 shows the cumulative LFs of Lya of the LAE 
candidates in the SSA22 field and in the SDF. There is a 
clear difference. As summarised in Table 4, the number den- 
sity of the z ~ 6.5 LAE candidates with L^ya > 3 x 10 42 
erg s _1 in the SSA22 is a factor of 2.5 smaller than that 
in the SDF. The number ratio (SDF/SSA22) is 3.5 for 
L Lya > 5 x 10 42 erg s _1 and 3.4 for L Lya > 7 x 10 42 erg 
s~ . In Figure 6, the SSA22 LF seems to be reproduced if 
the number density of the SDF LF is reduced by a factor 
of 0.3. On the other hand, a factor of 0.6 reduction of I/Lya 
in the SDF LF also provides an agreement with the SSA22 
LF. The reduction of LLya corresponds to a scenario with a 
different neutral fraction of hydrogen in the two fields. We 
will discuss this point more in section 5.2. 
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Figure 6. Cumulative Lyo luminosity functions (LFs) of z > 6 
LAEs. The thick solid and dot-dashed histograms are the LFs 
of photometric LAE candidates at z = 6.5 in the SSA22 field 
and in the SDF, respectively. The thin solid histogram is the LF 
of spectroscopic LAE in the SSA22 field. No correction of the 
completeness has been applied. The dotted histogram is the SDF 
LF whose number density is reduced by a factor of 0.3. The dashed 
histogram is the SDF LF with a reduction factor of 0.6 for the 
Lya luminosity. The thin solid curve is the LF for z = 5.7 LAEs 
in the SDF by Shimasaku et al. (2006). The open circle with an 
error-bar is the estimation based on one spectroscopic LAE at 
z = 7 in the SDF by Ota et al. (2008). 

4.2 Number density of i'-drop galaxies 

Here we compare the number densities of i'-drop galaxies in 
the SSA22 and in the SDF. This comparison may be use- 
ful to resolve the reason why the number of z ~ 6.5 LAE 
candidates in the SSA22 is so small relative to that in the 
SDF. Namely, the small number of the LAE candidates in 
the SSA22 may suggest that the field is just a void-like re- 
gion at z ~ 6.5. If so, the number density of i'-drop galaxies 
may be also small. However, we should be care about the 
difference of the depths along the line of sight for sampling 
i'-drop galaxies and the LAE candidates: ~ 250 Hq \ Mpc 
for i'-drop galaxies but 42 \ Mpc for the LAE candidates. 

To select i'-drop galaxies among objects detected in 
the z' band, we adopt the following criteria for the SSA22 
(SDF): 

(a) 25.50 < z' < 26.00 (5<r), 

(b) i' - z' > 1.5 or i' > 27.53 (27.84) (2a), 

(c) B > 27.89 (28.38), V > 27.66 (27.69), and R > 27.03 (27.76) 
(3a), 

where a indicates 1-er uncertainty of photometry in each 
bands and the values in parenth eses for the S DF. These 
criteria are sim ilar to those in Na gao et al.l (|2004r ) and 
lOta et all (|2005l ). We do not have near-infrared data and 
cannot adopt any res t-fram e UV colour crite r ia as done by 
iBouwens et"afl (|2003D and IShimasaku et all (|2005l ). Thus, 
our sample of the i'-drop objects is contaminated by Galac- 
tic M/L/T dwarfs . However, the contamination is < 10% 
for z > 25.5 AB (|Ota et al.ll2005h . The survey area in the 
SSA22 is the same as the LAE candidates. 

The resultant numbers of the i'-drop galaxies are 25 in 



the SSA22 and 36 in the SDF. The surface densities are 
0.047 ± 0.009 arcmin" 2 in the SSA22 and 0.036 ± 0.006 
arcmin -2 in the S DF both which ar e consistent with that in 
the SXDS field bv lOta et all |2005l ). Therefore, the number 
density of i'-drop galaxies in the SSA22 is 'normal'. This 
may suggest that a large difference of the number densities 
of the LAE candidates in the two fields is not caused by 
different number densities of dark matter haloes. However, 
we should caution ourselves that the sampling redshift of 
the i'-drop selection is different from that of the LAE se- 
lection. On the other hand, the spatial distribution of the 
i'-drop galaxies in the SSA22 seems to be well correlated to 
that of the LAE candidates as shown in Figure 3. This may 
indicate a similar redshift range of the i'-drop galaxies to 
the LAE candidates, while we should wait for spectroscopic 
confirmations to reach the final conclusion. 



4.3 Effect of smaller aperture size 

Because of finer FWHMs in our final images, the aperture 
sizes in our photometry are smaller than those in the SDF. 
This difference may affect the estimation of Lycv luminosity 
and the resultant LF unless the Lya emissions are point- 
like in our image. Indeed, some LAE candidates seem to be 
spatially extended as found in Figure 4. To examine this 
issue, we have performed the LAE selection again in the im- 
ages convolved with a Gaussian kernel so as that FWHMs 
of point-like sources become l."0 which is the same in the 
SDF images. The selection criteria are the same as §3, but 
the limiting magnitudes are ~ 0.3 mag shallower. The pho- 
tometry was made within a diameter of twice of the FWHM 
(i.e., 2."0). After removing objects affected by spiders of 
bright stars, we have found 6 objects with NB912 < 25.76 
(5-u) and 4 objects with L\, yoi > 5 x 10 42 erg s _1 . This num- 
ber is consistent with 3 objects found in the finer FWHM 
images other than the objects SSA22-z6p5-13 and -14 whose 
photometry affected by a foreground object. Therefore, the 
smaller aperture sizes due to finer FWHMs do not affect our 
conclusion. 



5 DISCUSSIONS 
5.1 Cosmic variance 

We have found a large difference in the number densities 
of 2 ~ 6.5 LAE candidates in the SSA22 and in the SDF. 
Here we examine if this difference can be explained by vari- 
ance of the number density of dark matter haloes, i.e. cosmic 
variance. For this aim, we assume the SDF LF to be repre- 
sentative. 

First, we calculate the expected number of z ~ 6.5 
LAE candidates in the SSA22 from the observed number 
in the SDF. As shown in the third column of Table 5, the 
expected numbers, A?c X p , are a factor of about 3 larger than 
the observed numbers, N a \>s. Next, we estimate a fractional 
standard deviation by cosmic variance, a cv , according to 
ISomerville et al] (|2004l ). As summarised in the fourth col- 
umn of Table 5, the a cv values become smaller for less lu- 
minous LAE candidates because the bias parameters are 
smaller. The survey volume of our SSA22 survey is smaller 
than that of the SDF, and thus, the resultant <t cv values are 



Table 5. Number of z ~ 6.5 LAE candidates in the SSA22 survey 
expected from the SDF luminosity function and cosmic variance. 





N oh3 a 


/V b 

^*cxp 


0"cv C 


P« N oba ) d 


(10 42 erg s" 1 ) 






(36) 


(%) 


> 10 





3.1 


49 


12 


> 7 


2 


6.8 


46 


14 


> 5 


5 


17.4 


41 


6.5 


> 3 


11 


34.7 


36 


6.6 



a Observed number of the LAE candidates in the SSA22 survey. 
b Number of the LAE candidates found in the SSA22 survey ex- 
pected from the SDF result. 

c Fractional sta ndard deviation by cosm ic variance estimated by 
the method in Somcrvillc et al. | d2004h based on the observed 
number density in the SDF and the SSA22 survey volume. 
d Probability to have a number equal to or smaller than the ob- 
served number. This is estimated from the expected number by a 
Monte-Carlo simulation taking into account the Gaussian cosmic 
variance and Poisson error. 

larger than those estimated in iKashikawa et all for 
the SDF. Although we ado pted the lines for z — 6 in figure 
3 of lSomerville et all <|2004h which are the closest to z = 6.5, 
the uncertainty caused by the redshift difference would be 
small. 

In order to examine the statistical significance of the 
smallness of the observed numbers, we have performed a 
Monte-Carlo simulation as follows; Assuming the cosmic 
variance to be Gaussian, we first draw a random number 
from the Gaussian distribution with its mean of A^xp and 
the standard deviation of N ex _ p a cv . Then, we draw another 
random number from the Poisson distribution with its pa- 
rameter of the first random number and store the second 
random number as the "expected observed number", . 
We repeat this procedure 100,000 times for each Lya lumi- 
nosity bin. Finally, we count the realizations whose 
is equal to or smaller than the real observed number N a b B - 
The resultant probabilities are listed in the last column of 
Table 5. We find that the probability to have A^bs of LAE 
candidates with L^ ya > 5 — 3 x 10 42 erg s _1 as small as 
in the SSA22 field is less than 7%, while the probability 
is higher for more luminous LAEs due to a larger Poisson 
error. From the Monte-Carlo simulation, we may conclude 
that the SSA22 field is a rare void at z ~ 6.5. 

However, we find no significant difference of the number 
density of i'-drop galaxies in the SSA22 and the SDF. In- 
deed, the number of i'-drop galaxies in the SSA22 expected 
from that in the SDF is consistent as summarised in Table 6. 
Note that the cr cv for this case is smaller than those in the 
LAE case because the survey volume is larger. On the other 
hand, the survey depth along the line of sight for i'-drop 
galaxies is about 6 times larger than that for LAEs. This 
point may explain the consistency of the number density of 
i'-drop galaxies even if there is a void at z ~ 6.5 traced by 
LAEs. 
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Table 6. Same as Table 4 but for i'-drop galaxies. 



z' 


A ba A exp 




Pi> N ohs ) 


(AB) 




(96) 


(%) 


25.5-26.0 


25 19.5 


27 


23 



ber density of the observable LAEs decreases in the SSA22. 
As found in Figure 6, the cumulative LF in the SSA22 is 
almost reproduced if the Lya luminosit y in the SDF is re- 
duced by a factor of 0.6. According to ISantosI l|2004h . the 
Lya transmission through the IGM is roughly proportional 
to 1/xhI: where xm is the hydrogen neutral fraction in the 
IGM. Therefore, the factor of 0.6 reduction of Lya luminos- 
ity can be interpreted as a factor of about 2 larger xhi in 
the SSA22 than in the SDF . 

Kashikaw a et al.l (|2006h have already found a significant 
reduction of the Lya luminosity at z = 6.5 relative to z — 
5.7 in the SDF and argue that the IGM is more ne utral at 
z = 6 .5 than z = 5.7. Their estimation based on ISantos! 
(|2004 ) model is xm < 0.45 at z = 6.50 In the SSA22, there 
may be a further factor of ~ 2 reduction of Lya luminosity. 
This is translated into xm < 0.9 at z — 6.5 in the SSA22. 
Therefore, more than half hydrogen may be still neutral at 
z = 6.5 in the SSA22 field. 

The furth est known LAE i s IOK-1 at z = 6.96 found 
in the SDF by live et al l i|2006D . IOta et al.l (|2008h have con- 
strained ihi at 2 = 7 in the SDF based on this one object: 
xhi = 0.1-0.6. The Lya luminosity of this object is 1.1 x 10 43 
erg s _1 . On the other hand, we could not find any LAE can- 
didates with L Lya > 1 x 10 43 erg s" 1 at z = 6.5 in the SSA22. 
The number density of z = 7 LAE in the SDF is ~ 3 x 10" 6 
Mpc -3 for Z/Lya > 1 x 10 43 erg s _1 which may match with 
an extrapolation of the SSA22 LF shown in Figure 6. This 
may imply that the SSA22 at z = 6.5 is as neutral as the 
SDF at z = 7. 

If the z = 6.5 IGM in the SSA22 field is more neutral 
than that in the SDF and Lya luminosity in the SSA22 field 
is reduced, an effect may appear in the Lya EW distribu- 
tion. Figure 7 shows the observed Lya EW as a function of 
UV luminosity for z = 6.5 photometric LAE samples in the 
SSA22 field and in the SDF. Although we expect that the 
EWs of the SSA22 sample are systematically smaller than 
those of the SDF sample, we can not find such a trend. On 
the other hand, the UV luminosities in the SSA22 field may 
be systematically smaller than those in the SDF, which may 
indicate that the small number of the LAEs in the SSA22 
field is due to a cosmic variance. However, in any case, we 
can not extract any conclusion due to the small statistics of 
the SSA22 sample. A future work with a much larger sample 
is required. 



5.2 Difference of neutral fraction and implications 
for reionization 

As an alternative scenario, if neutral hydrogen remains more 
in the SSA22 than in the SDF at z = 6.5, Lya emission is 
more heavily obscured in the SSA22, and thus, the num- 



Kashikawa et al. I 1120061 ) obtained an upper limit on xjji- The 
Lya luminosity reduction (or IGM transmission of Lya) could be 
accounted for by almost ion ized IGM if there is no galactic wind 
(see Fig. 25 in lSantoslliooi) . 
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as thi < 0.45 (Kashika wa et all [200 (f). Therefore, we may 
obtain xm < 0.9 at z = 6.5 in the SS A22. This i s simil ar to 
that at z = 7 in the SDF reported bv lOta et all i|200ST ). 

In any case, the smallness of the number density of the 
z — 6.5 LAE candidates in the SSA22 clearly shows that 
there is a large fluctuation of the LAE n umber den s ity at 
this hi gh redshift as alread y suggested by iHu et alj (I2004I. 
l2005h ; iHu fc Cowid (|2006h : IHu et all l|2010l ). lOuchi et all 
(2008) also reports a factor of 5 variation in the number 
densities of z = 5.7 LAEs among 5 fields-of-view of the S- 
Cam. Therefore, we need a much wider survey area than the 
surveys reported here and previously in order to obtain a ro- 
bust mean LF of z = 6.5 LAEs and a robust mean neutral 
fraction in the IGM at the redshift. 



Figure 7. Observed Ly« equivalent width as a function of UV 
luminosity. The squares are the z = 6.5 LAE candidates in the 
SSA22 field and the crosses are those in the SDF. The arrows 
indicate upper or lower limits for Luv or EW, respectively. 



6 CONCLUSIONS 

This paper presents results of a deep survey of z = 6.5 LAE 
candidates in the SSA22 field based on the imaging data 
collected with Subaru/S-Cam. We have selected the LAE 
candidates by the standard technique based on a narrow- 
band excess by an emission line and drop-outs at shorter 
wavelengths. The adopted narrowband is NB912 which can 
capture Lya emission at z = 6.46-6.57. We have also se- 
lected i'-drop galaxies for a comparison. These galaxies exist 
at a similar high redshift to but more extended range than 
the LAE candidates. 

A summary of the observational findings is as follows; 
The number of the LAE candidates with LL ya > 3 x 10 42 erg 
s _1 in the observed area of ~ 530 arcmin 2 in the SSA22 field 
is only 14. The co-moving number density of the LAE candi- 
da tes is a factor of 3 smal ler t han that in the SDF rep orted 
by iTaniguchi et all (|2005l ) and iKashikawa et all (|2006h . On 
the other hand, the number densities of i'-drop galaxies in 
the two fields agree within the Poisson error. 

We have considered two possibilities accounting for the 
small number density of the LAE candidates in the SSA22: 
(1) fluctuation due to the large-scale structure (i.e. so-called 
cosmic variance) and (2) fluctuation of the neutral fraction 
of hydrogen in the IGM. In the cosmic variance scenario, 
the probability to have the small number density as in the 
SSA22 field is 7% if we assume the SDF to be the cosmic 
average. Therefore, the SSA22 field at z — 6.5 may be a 
rare void. However, there is no evidence of such a void in 
the i'-drop galaxies although their redshift range is different 
from the LAE candidates. 

If the neutral fraction in the IGM is higher in the SSA22 
than in the SDF, a part of Lya emission is more heavily 
obscured in the SSA22. Indeed, a factor of about 2 reduction 
of Lya luminosity in the SDF provides a good fit of the 
SSA22 Lya LF. This is interpreted as a factor of about 2 
higher neutral fraction in the SSA22 field than in the SDF. 
Even in the SDF, z = 6.5 Lya LF indicates a factor of 
2 reduction of Lya luminosity relative to z = 5.7 and this 
reduction puts a constraint on the neutral fraction at z — 6.5 
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